Controlling crystallization and its absence: proteins, colloids and patchy models 
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The ability to control the crystallization behaviour (including its absence) of particles, be they 
biomolecules such as globular proteins, inorganic colloids, nanoparticles, or metal atoms in an al- 
loy, is of both fundamental and technological importance. Much can be learnt from the exquisite 
control that biological systems exert over the behaviour of proteins, where protein crystallization 
and aggregation are generally suppressed, but where in particular instances complex crystalline as- 
semblies can be formed that have a functional purpose. We also explore the insights that can be 
obtained from computational modelling, focussing on the subtle interplay between the interparticle 
interactions, the preferred local order and the resulting crystallization kinetics. In particular, we 
highlight the role played by "frustration", where there is an incompatibility between the preferred 
local order and the global crystalline order, using examples from atomic glass formers and model 
anisotropic particles. 



I. INTRODUCTION 

Controlling crystallization is a subject of considerable 
importance both from a fundamental and an applied 
perspective. Chemical physicists want to understand 
how the interparticle interactions determine the most 
favoured crystal structure and the ease with which crys- 
tallization can occur. Biochemists want to know the best 
recipe to crystallize the protein in which they are in- 
terested, so that they can then determine its structure. 
Nanotechnologists want to know how to get colloids or 
nanoparticles to self-assemble into a given target struc- 
ture, such as a diamond lattice because of its potential 
importance in photonics. Metallurgists want to be able 
to predict which alloys will most readily avoid crystal- 
lization, and instead form a metallic glass. 1 

In this paper we want to take a theoretical and com- 
puter simulation perspective on the factors that control 
crystallization, and its absence, in these kinds of sys- 
tems. So far in the literature, isotropic models have 
been the starting point for much of the theoretical work 
in this area, and although there have been consider- 
able successes, this approach has its limits. For exam- 
ple, the experimentally observed "crystallization slot" 
for globular proteins, where the interprotein interactions, 
as measured by the second virial coefficient, are suffi- 
ciently attractive to encourage crystallization, but not 
so attractive as to lead to irreversible aggregation^ has 
been rationalised using isotropic model potentials £ An- 
other fruitful idea is that if a protein solution lies near 
to a metastable critical point associated with separa- 
tion into protein-rich and solvent-rich phases, the asso- 
ciated concentration fluctuations could enhance crystal 



nucleation^ However, proteins are not isotropic, but are 
anisotropic both in their shape and in their interactions. 
One reflection of the latter is that protein crystals typi- 
cally have much lower packing fractions 5 than the close- 
packed structures that are favoured by isotropic poten- 
tials. 

Similarly, for colloids and nanoparticles there are 
now a considerable array of different crystals that have 
been obtained just using isotropic particles, particu- 
larly through the use of binary mixtures. For ex- 
ample, just by varying the relative sizes of the two 
types of particles, surprisingly complex crystals can 
be formedi&i Very recently, the use of particles with 
charges of opposite signs has led to colloidal 8 -^ and 
nanoparticl o 10 ' 11 ! 12 analogues of ionic crystals, even rela- 
tively low-density structures such as zinc blende?. 12 How- 
ever, there are limits to the structures that can be 
assembled from isotropic particles^ and so there has 
been considerable recent interest in developing meth- 
ods to produce colloidal particles that are anisotropic 
in shap o 14 ' 15 ! 16 or in their interactions i 15 i 17 ' 18 i 19 i 20 i 21 For 
all these reasons, there has been much theoretical inter- 
est in beginning to explore simple anisotropic models and 
their effects on crystallizatio n 22 i 23 i 24 ' 25 i 26 ' 27 i 28 ' 29 and self- 
assembly ^i 31 ^^ 

Here we provide a particular viewpoint on this broad 
area of controlling crystallization and its absence. In Sec- 
tion [TT] we outline some of our perspectives on protein 
crystallization, particularly why it is important to take 
into account the evolutionary origins of proteins when 
trying to understand their crystallization behaviour. In 
Section [TH] we explore how the concept of frustration can 
provide physical insight into the structural origins of the 
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glass-forming ability for two of the most common models 
for simulating supercooled liquids. Finally, in Scction HVl 
we introduce some of our recent results that use model 
patchy particles to explore how the geometry of the in- 
terparticle interactions can both influence and be used to 
control the crystallization and self-assembly behaviour. 



II. PROTEIN CRYSTALLIZATION 

The question concerning protein crystallization that 
we have been particularly seeking to address 3 ^ is: Why 
are globular proteins seemingly so hard to crystallize? 
Although it is not true of all proteins — lysozyme and 
insulin provide notable examples of proteins that crys- 
tallize relatively easily — there is plenty of anecdotal evi- 
dence that proteins are often very difficult to crystallize. 
Furthermore, the rise of structural genomics initiatives, 
has now made it possible to start to quantify this diffi- 
culty. For example, the success rate for producing X-ray 
quality crystals has been estimated to be roughly 20% for 
those prokaryotic proteins that can be expressed in sol- 
uble forni i 35 i 36 Moreover, the crystallization approaches 
used are mainly empirical. For example, screening kits 
allow large numbers of solution conditions, which have 
previously been found to be useful for crystallization, to 
be quickly tested. 

Our perspective on this question 3 ^ is that the need for 
proteins to function within the crowded cellular environ- 
ment places major constraints on the surface properties 
of proteins. In particular, not only must a protein inter- 
act correctly with its binding partners, but it must also 
make sure that it does not stick to anything else in the 
cell, including other copies of itself. Failure to do so is 
likely to be deleterious to the cell, as in the many protein 
aggregation diseases, of which the most relevant to this 
paper are those associated with native state aggregation 
and crystallization. For example, sickle cell anaemia in- 
volves the ordered aggregation of hemoglobin into fibrils, 
and some forms of anaemia and cataracts are the result of 
the crystallization of mutant forms of the haemoglobin 37 
and 7-crystallir* 3 ^ proteins, respectively. 

Thus, our hypothesis is that the surface properties of 
proteins have been selected in order to prevent native- 
state aggregation and crystallization in vivo (we term this 
'engineering out' of unwanted properties negative design) 
and that, because of the robustness of the mechanisms 
used, many proteins are difficult to crystallize even in 
the far-from-physiological conditions explored by the pro- 
tein crystallographer. Furthermore, there is a relatively 
simple way to test this hypothesis. If negative design 
is present, then random mutagenesis of surface amino 
acids should on average lead to a protein that is easier to 
crystallize, and the two such mutagenesis studies in the 
literature do indeed find such a correlationj 39 i 40 

More useful, though, would be to identify the mech- 
anisms used to achieve this negative design, since then 
this raises the possibility of developing rational strate- 



gies to overcome the negative design in order to make 
protein crystallization easier. Interestingly, bioinformatic 
analyses have provided a 'smoking gun' for the poten- 
tial role of surface lysine residues. Lysine is the most 
common surface amino acid, but yet is the most under- 
represented at functional interface a 41 i 42 and at contacts 
between proteins in crystals. 43 This of course begs the 
question: What is it doing there, if it is only reluctantly 
involved in interactions? One possible answer is that ly- 
sine helps to prevent unwanted interactions. 

Consistent with this suggestion, the group of Zygmunt 
Derewenda have shown that systematically mutating sur- 
face lysine residues, particularly to alanine^ almost in- 
variably enhances the crystallizability of a proteini 36 i 45 
The long lysine side-chain has substantial conformational 
entropy, which would be lost if constrained at a protein- 
protein interface. Thus, it has been suggested that ly- 
sine stabilizes the free protein in solution by providing 
a 'surface entropy shield'^ and replacing lysines with 
the compact alanine is expected to make crystallization 
thermodynamically more favourable. 

One possible objection to our negative design hypoth- 
esis is that, given that proteins are irregular objects sub- 
ject to dynamical fluctuations, it is not surprising that 
they are hard to crystallize. However, further evidence 
that the low crystallizability of most proteins is not an 
intrinsic property, but one that arises due to selection, 
comes from instances when protein crystallization does 
occur in vivo, because it achieves some functional pur- 
pose. Nature has no problem getting proteins to crys- 
tallize, even in the 'dirty' environment of the cell, when 
there is a positive selection pressure for this. 

There are many fascinating examples of this in vivo 
protein crystallization, as catalogued in a recent review 4£ 
Such crystallization represents a rather neglected area 
of biological self-assembly, and one which is worthy of 
further study. Here we give just a few examples. Certain 
genera of viruses that infect insect larvae, coopt the cell 
to express large quantities of specific proteins during the 
late stages of infection. These then crystallize around 
the viral particles to provide a protective environment 
for the viruses after the death of the insect larvae, as 
illustrated in Fig. |ljb) . On ingestion of the crystals by 
a new insect larva, the protein crystals dissolve in the 
alkaline environment of the gut, releasing the virus to 
infect the new host4£ 

A number of examples of crystalline proteins in the cell 
correspond to proteins stored in secretory granules. Fig. 
[TJa) illustrates one particularly specialized example from 
the protist (a type of single-celled eukaryotic organism) 
Paramecium^ Large numbers of these granules, which 
are called trichocysts, are attached to the outer mem- 
brane of this organism, and are assembled from three 
different families of closely related polypeptides, each lo- 
calized to different regions of the trichocyst to form a 
structure that is of the order of 3-4 /im in size. In re- 
sponse to an external stimulus the trichocysts can be 
simultaneously extruded from the cell. In this process, 
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FIG. 1: Some examples of protein crystallization occurring in 
the cell, (a) A trichocyst attached to the outer membrane of 
Paramecium, (b) An encapsulated rod virus of the granulosis 
virus of Plodia interpunctella. (c) Paracrystalline arrays of 
an iridescent virus in the epidermis of a larva of Culicoides 
variipennis sonorensis. The scale bars in (b) and (c) corre- 
spond to 0.1 /im and 1 M m i respectively. Reproduced with 
permission from Refs. (a) [4g and (b) [4S| and (c) |50| . 



the crystalline trichocysts expand in length by a factor of 
roughly eight. The purpose seems to be defensive. The 
explosive release of the trichocysts can push Paramecium 
away from a potential predator, giving it a chance to 
escaped 

In both these examples, a high level of control is as- 
serted over the crystallization process in order to deter- 
mine the location, size and shape of the crystal. The 
mechanical properties are also important for their func- 
tion with the virus-encapsulating crystals required to be 
tough in order to provide a protective environment, and 
the trichocysts able to undergo an irreversible expansion 
triggered by the presence of calcium ions. 

Our final example is of the crystallization of virus par- 
ticles themselves (Fig. [He)), and has perhaps more to 
do with an absence of a selection pressure to prevent 
crystallization (why should it bother the virus?) than a 
functional purpose for the crystals. The most dramatic 
examples come from iridoviruses, which have large icosa- 
hedral capsids. This type of virus was first identified due 
to the iridescent colours imparted to insect larvae (hence 
the name)^ due to the Bragg scattering of visible light 
by the crystals^ 



III. GLASS FORMATION 

There has been much work trying to understand the 
unusual dynamic properties of supercooled liquids as the 



glass transition approached^ To do this from a simu- 
lation perspective, model systems have been developed 
that are robust glass-formers and show no tendency to 
crystallize. Such models also allow one to probe the 
structural determinants of a system's glass-forming abil- 
ity. A key notion is that of 'frustration' ) 55 ^ 6 which is 
said to occur when the preferred local order is incompat- 
ible with global crystalline order. Typically, the liquid 
structure reflects the preferred local order, and so when 
frustration is present, significant structural reorganiza- 
tion will be necessary to nucleate a crystalline phase. 

To apply these ideas, one needs to be able to identify 
the preferred local order, but, because of the inherent 
disorder associated with a liquid, sometimes this can be 
hard to achieve directly from analysing liquid configura- 
tions. An alternative approach is to instead look at the 
structures of isolated clusters, since such a cluster can 
adopt the preferred structure without having to mould it- 
self to any surrounding medium. This approach was first 
applied by Charles Frank to rationalize why small liquid 
metal droplets could be substantially supercooled^ 

One of the potentials most commonly used by simula- 
tors is the Lennard- Jones potential. However, the one- 
component Lennard- Jones fluid is not a good system to 
look at the behaviour of supercooled liquids, since crys- 
tallization into a close-packed solid occurs relatively eas- 
ily. To generate a good glass-former, Dzugutov modi- 
fied this potential by adding a barrier in the potential at 
roughly y2 times the position of the potential minimum 
(Fig. UK a)), hence energetically disfavouring close-packed 
crystals because of their octahedral interstices^ Instead, 
polytetrahedral order^ and local icosahedral order in 
particular, is preferred. As can be seen in Fig. UK a), the 
distances present in the 13-atom icosahedral cluster avoid 
sampling this bump in the potential. 

Even though there is a stable body-centred-cubic 
phase,— the resulting Dzugutov liquid can be easily su- 
percooled. Interestingly, the Dzugutov system is one of 
the only liquids in which an increase in a structural length 
scale has been detected on increased supercooling, in par- 
ticular domains with local icosahedral order grow in size 
as the temperature is decreased^ However, these do- 
mains are not compact, but are ramified in structure. 
These results raise the questions: What are the origins of 
the non-compactness, and might this shape be a generic 
property of domains of enhanced local order within su- 
percooled liquids? An analysis of the structures of iso- 
lated clusters is able to provide some answers. This type 
of ordering is also clearly seen in the clusters, but an anal- 
ysis of the energetics shows that it is a consequence of the 
unusual shape of the Dzugutov potential. Therefore, this 
domain shape is unlikely to be universal. 

The lowest-energy Dzugutov clusters are all aggre- 
gates of interpenetrating and face-sharing icosahedra 
(Fig. HKb)).— As the cluster size increases first chains, 
then discs, then rings, and finally 3-dimensional porous 
networks of icosahedra are seen. The energetic causes for 
this ordering are also clear. The bump in the potential 
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FIG. 2: (a) The Dzugutov (Dz) and Lennard- Jones (LJ) po- 
tentials compared. The vertical lines show the positions of the 
pair distances in the 13-atom icosahedron and their heights 
are proportional to the numbers of pairs with that distance, 
(b) Some of the particularly stable structures for Dzugutov 
clusters. Each cluster is labelled by its size and point group. 




promotes local icosahedral order. However, icosahedral 
structures are inherently strained, e.g. the distance be- 
tween adjacent atoms on the surface of the regurlar 13- 
atom icosahedron is 5% longer than that to the central 
atom (Fig. Ufa)). This strain, and the associated ener- 
getic penalty, grows rapidly for compact polyicosahcdral 
structures. The resulting non-compact icosahedral ag- 
gregates represent a compromise that maintains the lo- 
cal icosahedral coordination, but which does not involve 
excessive strains. Similar behaviour has also been seen 
for variants of the Dzugutov potential^ 

A second approach used to generate a good glass- 
former is to introduce two atom types. Partly this is 
because the formation of a compositionally ordered crys- 
tal can require the atoms to diffuse significantly further to 
find the correct environment than for the one-component 
case. However, there is more to it than this, since the in- 
teractions between the particles and the composition also 
need to be tuned to reach those regions of the parameter 
space where crystallization is particularly difficult i 64 i 65 
The most commonly used model of this type is the Kob- 
Andersen binary Lennard-Jones mixture^ i.e. the po- 
tential is 



VBLj(Vy) = 4e Q/3 
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(3.1) 



where a and j3 are the atom types of atoms i and j, 
and e a /3 and 2 1 / e a a p are the pair well depth and equi- 
librium pair separation, respectively, for the interaction 
between atoms i and j. The Kob-Andersen parameters, 



FIG. 3: (a) The energies of the global minima for four series 
of Kob-Andersen BLJ clusters with different numbers of B 
atoms. To make particularly stable clusters stand out, the 
energies are measured with respect to E ave , a fit to the en- 
ergies of the clusters. For clarity, the curves have also be 
displaced with respect to each other, (b) Some of the par- 
ticularly stable structures identified in (a) with all B atoms 
completely coordinated, (c) A particularly stable cluster at 
near equimolar composition. 



a A B = 0.8 <taa, &BB = 0.88 cr A a, <^AB = 1-5 c AA , and 
(bb = 0.5 e A A, are non- additive and strongly favour mix- 
ing. At the canonical composition A4B the ground state 
is a coexisting pure A face-centred-cubic crystal, and an 
AB CsCl-type crystal^ with low-lying layered A4B crys- 
tals also possible^ However, crystallization has never 
been seen in a simulation. 

Here, we use isolated clusters to analyse the preferred 
coordination environments of A atoms around B atoms, 
since at the glass-forming composition B atoms are in 
a minority. In particular, we locate the particular sta- 
ble clusters with one, two, three and four B atoms to 
see how these preferred environments can organize into 
larger structures (Fig. [3]). For coordination of a single B 
atom, clusters with 8, 9 and 10 atoms have similar stabil- 
ity. All three structures are based on a square antiprism 
of A atoms, but with possible capping atoms over the 
two square faces. This stability of the square antiprism 
is also seen for AgB2 and AgB3, where the additional 
capping atoms are now B atoms. One way of combining 
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such coordination environments is illustrated by A12B2 
and A 16 B 3 and involves the sharing of the square-faces 
of the square antiprisms to form linear aggregates. Simi- 
lar columns of square anti-prisms are found in the AI2CU 
crystal. A18B3 illustrates another possibility and involves 
the sharing of two triangular faces of A12B2 with a bi- 
capped square antiprism. 

Another possible coordination environment for the B 
atoms is the tricapped trigonal prism. Although this 
structure is not stable for A9B (it can relax to a mono- 
capped square antiprism by a single diamond-square 
process 6 -^) it becomes a common environment for larger 
clusters. This is illustrated by the structures A15B3, 
A 17 B 4 and A 2 oB 4 , which have one, two and four tri- 
capped trigonal prismatic environments, respectively, al- 
beit with B atoms in some of the capping sites. 

These results tie in well with studies of the local 
structure in liquid configurations, which have found that 
square antiprism, and trigonal prismatic environments 
increasingly dominate the local structure in the liquid as 
the temperature is decreased^ 

So how do these results help to rationalize the sys- 
tem's ability to avoid crystallization. Firstly, these pre- 
ferred local environments are not present in the lowest- 
energy crystal structures. However, there are still a 
number of crystalline structures that are significantly 
lower in energy than the liquid that do involve these 
environments. 64 ' 65 ' 69 Secondly, the diversity of environ- 
ments and ways that these can pack is likely to frustrate 
the formation of a uniform crystal. 

In contrast to the A4B system, at equimolar compo- 
sitions the Kob- Andersen BLJ model easily crystallizes 
into a CsCl-type structured Again, isolated clusters can 
help to understand this behaviour. Examining AsB^ 
clusters, we found AgBy to be particularly stable and 
to exhibit the bulk crystal structure (Fig. [3th)), hence 
showing the absence of frustration at this composition. 




FIG. 4: The geometry of the interactions between the model 
particles. In this example, there are four patches regularly 
arranged on the discs, with their directions described by the 
patch vectors, pi. Patch 4 on particle i interacts with patch 
2 on particle j because they are closest to the interparticle 
vector. 



that is very much in its infancy with much to be discov- 
ered. Our interest is to probe how the anisotropy of the 
interparticle interactions can be used to control a sys- 
tem's crystallization and self-assembly behaviour, and, 
in particular, through its effect on the degree of frustra- 
tion. As well as being of fundamental interest, the hope 
is that the results will also be useful for understanding 
the crystallization behaviour of proteins and anisotropic 
colloids. 



A. Potential 

Here, we model the patchy particles with a single-site 
potential, i.e. each particle is represented by just a single 
site, but which has both position and orientation. The 
potential consists of an isotropic repulsion, which is based 
on the Lennard-Jones potential 



IV. PATCHY MODELS 



V LJ (r)=4e 



(4.2) 



In the above models with isotropic interactions, the 
connection between the preferred local structure and the 
form and parameters of the potential can be quite sub- 
tle. To probe the connection between crystallizability 
and local structure further, it would be useful to have 
a model where the local structure can be more directly 
controlled. To achieve this goal in a simple and flexible 
way, here we use 'patchy' particles, where the particles 
are spherical and only interact strongly when the patches 
on adjacent particles are aligned, but we should note that 
other ways of using anisotropic interactions to control the 
local structure have recently been explore d 70 ' 71 . There 
has been much recent interest in such patchy models from 
the perspective of protein crystallization , 22 ' 23 ^^ 6 -^ self- 
assembl y 29 ' 30 (particularly into monodisperse clusters in 
a way similar to virus self-assembl y 30 ' 31 ' 33 ) and the dy- 
namics of supercooled liquids ) 72 ' 73 but it is still an area 



but where the attraction is modulated by an orientational 
dependent term, Va llg . Thus, the complete potential is 



v-(r- n no - I Vu{r ^ r< 



r < <xlj 
(4.3) 

where Oj is the orientation of particle i. The patches are 
specified by a set of patch vectors, {pi}, as illustrated in 
Fig. [4] V" an g has the form. 




where a is the standard deviation of the Gaussian, akij 
is the angle between the patch vector k on particle i and 
the interparticle vector , and /c m in is the patch that 
minimizes the magnitude of this angle. Hence, only the 
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patches on each particle that are closest to the interpar- 
ticle axis interact with each other, and the potential is 
continuous as a function of the orientations of the parti- 
cles. V ang — 1 when two patches point directly at each 
other, but falls off as the patches deviate further from 
the perfect alignment. 

One of the nice features of this single-site potential 
is that, given a set of patch vectors, the potential has 
only one parameter er, which determines the widths of 
the patches. Furthermore, as l/a ^ the isotropic 
Lennard-Jones potential is recovered. Hence, it is pos- 
sible to systematically study the behaviour of the model 
as a function of the degree of anisotropy with the well- 
characterized Lennard-Jones model as one limit. 



B. Two-dimensional crystallization 

We first discuss the application of this model to crys- 
tallization in two dimensions to illustrate the effects of 
the geometrical arrangement of the patches, since visu- 
alization is easier than in three dimensions. We choose 
to study particles with 4, 5 or 6 patches arranged regu- 
larly on the surface of the disks. Particularly interesting 
will be the 5-patch system as the local five-fold symmetry 
of the patches is incommensurate with global crystalline 
order. 

The two-dimensional Lennard-Jones reference system 
crystallizes into a close-packed crystal with a hexagonal 
arrangement of neighbours around each particle. Unsur- 
prisingly, for particles with a regular hexagonal array of 
patches, the anisotropy reinforces this behaviour. Crys- 
tallization is easy (Fig. 02c)) and the close-packed crystal 
is lowest in energy for any combination of pressure and 
a. 

In the 4-patch system there is competition between a 
low-density square crystal in which each patch points di- 
rectly at a nearest neighbour, and a higher-density hexag- 
onal crystal. The latter is oricntationally ordered with 
two opposite patches pointing at nearest neighbours, and 
the other two at next neighbours. As the square crystal 
has more pairs of nearest-neighbour patches pointing di- 
rectly at each other, it is energetically preferred when 
the patches are sufficiently narrow (Fig. [5(a)). However, 
because of its more open structure, it becomes desta- 
bilised relative to the hexagonal crystal as the pressure 
is increased. Crystallization to the square crystal occurs 
relatively easily in the appropriate region of the phase 
diagram (Fig. [D(a)). 

For the 5-patch system the situation is more complex. 
No crystal phases are possible where all the patches point 
directly at neighbouring particles. For example, on cool- 
ing this system in the region of the phase diagram where 
the patches prefer to be aligned, the resulting configu- 
ration has no overall crystalline order. However, most 
of the particles have one of the two local environments 
highlighted in Fig. G3b). From these two motifs, two 
crystals can be constructed where every particle has an 



identical environment (Fig. [6(b)). In both crystals each 
particle has five nearest neighbours, and the crystals can 
be thought of in terms of tilings of squares and equilat- 
eral triangles, where every vertex is surrounded by three 
triangles and two squares. The nomenclature for these 
semi-regular tilings is (3 2 .4.3.4) and (3 3 .4 2 ), and derives 
from the ordering of polygons around each vertex^ 4 How- 
ever, we will use the notation 'er' and 'H' by analogy to 
the Frank-Kasper phases of these names that can be en- 
visaged as square-triangle tilings in two of their three 
dimensions^ 

When the patches are sufficiently narrow, the H- and 
(T-phases represent the best compromise between the five- 
fold symmetry of the particles and crystalline order. 
However, the patches cannot point directly at each other 
in these crystals, and so the system is frustrated. The 
mean deviation of the patches away from the nearest- 
neighbour interparticle vectors is smaller for the c-phase 
(7.2°) than for the H-phase (9.6°), making the former 
lowest in energy for intermediate values of a at low pres- 
sure. However, at small cr, the energetic penalty as- 
sociated with these deviations becomes significant, and 
causes the H-phase to distort so that three of the five 
patches can point directly at each other. This distortion 
reduces the aspect ratio of the unit cell and the triangles 
deviate significantly from the ideal equilateral geometry. 
No such similar distortion is possible for the er-phase, and 
so the H-phase becomes lowest in energy at small a (Fig. 
Mb))- 

The H- and cr-phase represent only two of the possible 
tilings of squares and triangles. These tilings need not 
be crystalline, and dodecagonal quasicrystalline pack- 
ings are also possible. Indeed, these quasicrystals have 
been seen in both metallic alloys^ and macromolecu- 
lar systems^ under circumstances near to where crys- 
talline square-triangle phases are stable. Therefore, a 
stable quasicrystalline state remains an intriguing possi- 
bility for the current system. 

The more dense hexagonal crystal is most stable at 
higher pressure, and closer to the isotropic limit. Al- 
though positionally ordered, the crystal is orientationally 
disordered, because of the incompatibility of the five-fold 
symmetry of the particles and the six-fold symmetry of 
the crystal (Fig. Mb))- 

An interesting comparison to the current system, 
which has five-fold symmetry in the attractive interac- 
tions, is provided by a study of hard pentagons^ Al- 
though the shape of these particles does frustrate crystal- 
lization to some extent, the lattice formed by the particle 
centres in the only stable orientationally-ordered crys- 
talline phase, is based on a slightly distorted hexagonal 
lattice. 

Also relevant is a recent study of a two-dimensional 
system of particles interacting with a Lennard-Jones po- 
tential plus an anisotropic term that favours the forma- 
tion of isolated five-fold rings of particles^. The presence 
of such order again frustrates crystallization, and, inter- 
estingly, on varying the strength of the anisotropic term 
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FIG. 5: Structures resulting from Monte Carlo cooling simulations for two-dimensional particles with (a) four (b) five and (c) 
six regularly spaced patches. In each case, there are 100 particles and a — n/12 ~ 0.262. The pressures are (a) 0.1, (b) 0.2 
and (c) 1.0e/<7Lj 3 . The dashed lines show the periodic boundary conditions. In (a) and (c) crystallization clearly occurs. The 
crystals are not quite perfect, because unless crystallization occurs at the correct orientation with respect to the boundary 
conditions and with the correct number of lattice points in each direction, there will not be the right number of particles to 
form a perfect crystal, and so some defects are necessarily present. In (b) there is no overall crystalline order, but two common 
motifs are highlighted. 



in the potential, the dynamic properties of the liquid, 
such as the fragility, can be changed substantially. 



C. Three-dimensional crystallization 

In our applications of this model to three dimensions, 
we again look at how low-density crystalline phases can 
be stabilized by the patchy interactions. The two sys- 
tems that we consider are 4- and 6-patch particles with 
a regular tetrahedral and octahedral arrangement of the 
patches, respectively. For these systems, the stable crys- 
talline phases at sufficiently low pressure and a would be 
expected to be a diamond and a simple cubic lattice, re- 
spectively, with each patch directly point at a neighbour- 
ing particle. These have maximum packing fractions of 
34% and 52% compared to 74% for the face-centred-cubic 
crystal favoured by the isotropic Lennard- Jones poten- 
tial. 

One might have expected these two systems to behave 
quite similarly in the regime where the open structures 
are preferred, but in fact their crystallization behaviour 
is quite different. The 6-patch system is able to crystal- 
lize easily, leading to a step-like decrease in the energy 
on crystallization (Fig. Ufa)), sometimes even giving a 
perfect defect-free crystal. By contrast for the 4-patch 
system at best only partial crystallization is observed, 
and on ordering the energy does not exhibit a step, but 
instead changes continuously. 

As in Section IIIH examining the structure of isolated 
clusters can give us some clues to the physical origins 
of this behaviour, because the clusters provide a picture 
of the preferred local order for the system. For both 



systems, when the patches are sufficiently narrow, as 
expected, the lowest-energy clusters exhibit open struc- 
tures where the maximum coordination number is equal 
to the number of patches. For the octahedral particles, 
for a < 0.5 most of the clusters are cuboidal nanocrys- 
tals with the simple cubic structure (e.g. the 3x3x3 
cube in Fig. 0c)). The system is unfrustrated and it is 
unsurprising that crystallization is relatively easy. 

By contrast, the global minima of the tetrahedral sys- 
tem for 0.2 < a < 0.45 do not exhibit the structure of 
the stable crystalline phase, but instead are based on do- 
decahedral cages (the 20-particle dodecahedron is shown 
in Fig. Efb)). Only for a < 0.2 are clusters with the 
diamond structures lowest in energy, e.g. the 26-particle 
structure in Fig. EJb) . In diamond the particles form six- 
fold rings, but which are puckered, whereas the dodeca- 
hedra are characterized by planar five-fold rings, where 
the 108° bond angles are very close to the 109.57° angles 
between the tetrahedral patches. The effect of this differ- 
ence is small until the patches become very narrow, and 
so for intermediate values of a, the dodecahedral clus- 
ters are preferred because they have a greater number of 
bonds, e.g. for the 20-atom dodecahedron there is only 
one unused patch per particle, whereas the particles on 
the surface of diamond clusters often have two unused 
patches. Thus, the results for these clusters suggest that 
the tendency of the liquid to form structures involving 
five-fold rings is one of the reason underlying the much 
greater difficulty of crystallization in the tetrahedral sys- 
tem. 

Furthermore, as well as this competition between five- 
and six-fold rings, in this model the diamond and hexago- 
nal diamond (the one-component equivalent of wurtzite) 
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FIG. 6: Structural phase diagram showing how the ground- 
state structure in two dimensions depends on patch width and 
pressure (in units of e/oLj 3 ) for particles with (a) four and (b) 
five regularly spaced patches. The unit cells of the crystalline 
structures are depicted with thick lines, and the dashed lines 
show the periodic boundary conditions. 



structures are practically degenerate. The structural dif- 
ference between these two crystals is that all the six-fold 
rings in diamond have a form analogous to the chair iso- 
mer of cyclohexane, whereas in hexagonal diamond, some 
of these rings are analogous to the boat form of cyclohex- 
ane. Again, this further variety of local structural forms 
is not going to aid crystallization. 

These results are consistent with the work of Zhang 
et al. who attempted to crystallize diamond using sim- 
ilar patchy particles with the same tetrahedral patch 
geometry^ Crystallization to the diamond structure 
only reliably occurred when a seed crystal was intro- 
duced, or when an additional term in the potential 
was added that favoured staggered over eclipsed tor- 
sional configurations, and hence disfavoured hexagonal 
diamond. 

This work on model anisotropic particles is particu- 
larly timely given the rapid recent advances in synthc- 
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FIG. 7: (a) The caloric curves for cooling the 4-patch tetrahe- 
dral particles and the 6-patch octahedral particles from high 
temperature. The number of particles is 512, the patch width 
a is 0.3 radians, and the pressure is 0.1 e/aLj 3 . The resulting 
simple cubic crystal for the 6-patch system is also illustrated 
(the dots represent the corners of the simulation box), (b) 
and (c) Some of the global minima of isolated clusters for the 
4- and 6-patch systems, respectively. Each cluster is labelled 
by the number of atoms and point group. 



sizing anisotropic colloids ; 14 ' 15 ! 16 and because one target 
for research in this field is to get particles with tetrahe- 
dral symmetry to assemble into a diamond-like crystal^ 
Our results indicate that the crystallization of such tetra- 
hedral colloids might not be so straightforward, because 
of the potentially frustrating effects of the variety of lo- 
cal structures possible in the liquid phase. However, it 
is not clear how an additional potential term similar to 
that used by Zhang et al. could be introduced into the 
intercolloidal potential to alleviate this. Instead, an al- 
ternative strategy would be to use a binary system of 
oppositely-charged tetrahedral colloids, as this would pe- 
nalize the formation of rings with an odd number of par- 
ticles. In particular, this change would reduce the ten- 
dency for pentagons to form, hence removing some of the 
frustration and making crystallization more likely. 

An interesting comparison to the present results are 
provided by a recent study that took the Stillinger- 
Weber silicon potential 7 ^ and varied the strength of the 
anisotropic 3-body term, finding that it had a significant 
effect on the system's glass-forming ability 71 . For the 
original silicon potential crystallization into a diamond 
structure occurred relatively easily, but as the strength 
of the 3-body term is decreased the system becomes a 



9 



good glass-former in the region of parameter space where 
there is a crossover in the stability of the diamond and 
body-centred-cubic crystals, and it has been suggested 
that this is partially due to the structural dissimilarity 
between the liquid and the possible crystals. Similarly, 
it would be interesting to look at how varying the patch 
width affected the dynamics of crystallization for our sys- 
tems. 



V. CONCLUSIONS 

Through the examples of crystallization (and its ab- 
sence) considered in this paper, we hope to have shown 
the subtle interplay between the interparticle potential, 
the preferred local structure and the kinetics of crystal- 
lization. In particular, if the interactions favour a local 
structure that is incompatible with the global crystalline 
order, or almost equally favour a variety of different lo- 
cal environments, then crystallization is likely to be frus- 
trated. Furthermore, we have provided further empirical 
support that examining the structures of isolated clus- 
ters, as first done by Frank, 57 can potentially provide a 
clearer picture of the preferred local order, and hence the 
possible presence of frustration. 

An important question concerning the type of struc- 
tural approach advocated here is how it relates to other 
approaches to understanding the kinetics of crystalliza- 
tion, such as classical nucleation theory. In the latter, 
one of the key parameters in determining the ease of nu- 
cleation is the surface free energy associated with the 
crystal-liquid interface. The connection to the current 
approach is that this interfacial free energy is particu- 
larly sensitive to the degree of structural dissimilarity 
between the crystal and the liquid, i.e. the greater the 
structurally dissimilarity, the greater the interfacial free 
energy, and hence the harder crystallization becomes. 

In such instances where the liquid and crystal are struc- 
turally very different, one way of circumventing the large 
free-energy barriers to direct nucleation of the crystal is 
to undergo a two-step nucleation process, i.e. first nucle- 
ate a metastable phase, and then nucleate the final crys- 
tal form from within that phas o 80 ' 81 . For example, for 
the Stillinger- Weber silicon potential mentioned earlier, 
crystallization to the diamond structure only occurs at 
temperatures below that for a liquid-liquid phase tran- 
sition, which leads to the nucleation of a lower-density 
liquid phase that more closely resembles the solidil. 

The patchy models introduced in Section[lV]are partic- 
ular useful for studying crystallization, since the geome- 
try of the patches allows the system's structural propen- 
sities to be directly controlled in a simple and flexible 



manner. One of our original intentions for these models 
was to also use them to illuminate aspects of the crys- 
tallization of proteins. Therefore, a key question is how 
relevant these models are to real proteins. Clearly, the 
interactions between proteins are strongly anisotropic in 
character, i.e. in order for two proteins (in their native 
state) to come together, the regions of their respective 
surfaces that come into contact, the "patches" if you 
like, must be correctly aligned and oriented. Further- 
more, the patchy models, like proteins, are also able to 
form relatively low-density crystals. However, actual in- 
terprotein interactions are much more complex — proteins 
are anisotropic in shape, their surfaces are very hetero- 
geneous, and the interactions can depend sensitively on 
solution conditions, such as pH and the concentrations 
of other ions. 82 Thus, how much of this complexity one 
needs to capture in order to reproduce or gain insight into 
the crystallization of proteins is an open, but important, 
question. 

The examples described in Section |TT] illustrate the 
exquisite control that biological systems are able to exert 
over the interaction properties of proteins, in particular 
their ability to avoid native state aggregation or crystal- 
lization in the dense cellular environment, and the ability 
of particular proteins to form into complex crystalline as- 
semblies. There are some indications of how this control 
is achieved, e.g. the role of lysine residues in prevent- 
ing unwanted interactions, but there is much still to be 
learnt. Such knowledge will be particularly important if 
rational methods to make protein crystallization in vitro 
easier are to be developed that aim to overcome the evo- 
lutionary selection of the surface properties of proteins 
to prevent native state aggregation. 

In some ways, our patchy particles are probably some- 
what closer to some of the anisotropic colloids now be- 
ginning to be generated ] 14 ' 15 ? 16 Therefore, it will be par- 
ticularly interesting, once it becomes possible to generate 
such colloids in sufficient number and quality, and with 
their surfaces appropriately functionalized, to be able to 
probe their phase behaviour and phase transformation 
kinetics. Hopefully, some of the insights gained from our 
patchy models will help to guide the experimentalists in 
their quest to create crystals with useful photonic prop- 
erties. 
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